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ABSTRACT 

To investigate the nature of starbursts' dust, we constructed a model of the stars 
and dust in starburst galaxies and applied it to 30 observed starburst spectral energy 
distributions (SEDs). The starburst model was constructed by combining two stellar 
evolutionary synthesis models with a model describing the radiative transfer of stellar 
photons through dust. The stellar evolutionary synthesis models were used to compute 
the dust-free SEDs for stellar populations with ages between 1 x 10 6 and 15 x 10 9 years. 
Using a Monte Carlo radiative transfer model, the effects of dust were computed for aver- 
age Milky Way (MW) and Small Magellanic Cloud (SMC) dust, two different star/dust 
geometries, and locally homogeneous or clumpy dust. Using color-color plots, the star- 
burst model was used to interpret the behavior of 30 starbursts with aperture-matched 
UV and optical SEDs (and IR for 19 of the 30) from previous studies. From the color- 
color plots, it was evident that the dust in starbursts has an extinction curve lacking a 
2175 A bump, like the SMC curve, and a steep far-UV rise, intermediate between the 
MW and SMC curves. The star /dust geometry which is able to explain the distribution 
of the 30 starbursts in various color-color plots has an inner dust-free sphere of stars 
surrounded by an outer star- free shell of clumpy dust. When combined with other work 
from the literature on the Orion region and the 30 Dor region of the Large Magellanic 
Cloud, this work implies a trend in dust properties with star formation intensity. 

Subject headings: dust, extinction — galaxies: ISM — galaxies: starburst — galaxies: 
stellar content — radiative transfer 



1. Introduction 

The intrinsic spectral energy distribution (SED) of 
a galaxy can tell us much about that galaxy's stel- 
lar evolutionary history, but dust affects the intrinsic 
SED in a complicated fashion. The effects of dust 
on a galaxy's SED are determined by the physical 
properties and spatial distribution of the dust. Our 
knowledge of the physical properties and spatial dis- 
tribution of dust in galaxies other than the Milky 
Way (MW) and the Magellanic Clouds is very limited. 
Starburst galaxies are good objects in which to study 
dust as their energetic output is dominated by mas- 
sive stars and they are bright in the ultraviolet (UV) 
where the effects of dust are large. In addition, proper 
removal of the dust's effects is necessary in order to 
study the starburst's stellar evolutionary history (i.e. 
age, star formation rate, initial mass function, etc.). 

The physical properties of dust affecting a galaxy's 
SED are parameterized by the wavelength depen- 
dence of the dust's extinction, albedo, and scattering 
phase function. The extinction curve is the domi- 
nant physical property and the easiest to determine 
observationally. Extinction curves for stars in our 
Galaxy, the Small & Large Magellanic Clouds (SMC 
& LMC), and (tentatively) M31 have been deter- 
mined and found to have a common shape, but ex- 
hibit significant variations. Generally the extinction 
increases with decreasing wavelength, making the ef- 
fects of dust strongest in the UV. The most prominent 
feature in UV/optical extinction curves is the 2175 A 
bump, which varies greatly in strength in our Galaxy 
dWitt, Bohlin, fe Stecher 1984|; |Cardelli, Clayton, fe 



Mathis 1989| ) and the LMC jFitzpatrick 1985|), is ab- 



sen| b in the SMC ■{ Prcvot ct al. 198 4 , but sec Lcqucux 



d. 1982 for an exception), and seems weaker in 



M3 *l than in our Galaxy (Dianchi ct al. 1996). The 



behavior of the 2175 A bump in the LMC seems par- 
ticularly relevant to the question of the type of dust 
present in starbursts. The average extinction curve 
for the 30 Dor region shows a very weak 2175 A bump 
and a strong far-UV rise, while the average extinc- 
tion curve for stars outside the 30 Dor region shows a 
2175 A bump strength and far-UV rise comparable to 
that of the average Galactic 2175 A bump (playton 
& Martin 1985] |Fitzpatrick 1985| ; |Fitzpatrick 1986[ ). 
This behavior implies the extinction curve in star- 
bursts may be more like the extinction in the 30 Dor 
region, lacking a significant 2175 A bump and having 
a strong far-UV rise. In fact, the 30 Dor region has 



been described as the Rosetta Stone for starbursts as 
it is the most massive star cluster in the Local Group 
(Hunter et al. 1995) and it is close enough to study 
individual stars (Walborn 1991). 

The spatial distribution of the dust plays a cru- 
cial role in determining the strength of dust's effects 
on the SED of a galaxy. Witt, Thronson, & Ca- 



puano (1992, hereafter WTC) investigated the effects 
of different spherical distributions of dust and stars 
(termed galactic environments) using a Monte Carlo 
radiative transfer model. Their finding was that the 
SED of a galaxy will always be dominated by the least 
attenuated stars as they contribute the most flux to 
the SED. As a result, the signatures of dust seen in the 
galaxy SED will almost always imply small amounts 
of dust. Recently, Witt & Gordon (1996) studied the 
radiative transfer in a clumpy (two-phase) dust distri- 
bution. They found the dumpiness resulted in more 
photons escaping than in the equivalent homogeneous 
case as the inter-clump medium provides small opti- 
cal depth paths. The dust in galaxies is distributed 
in clouds (e.g. [Dickey fc Garwood 1989| ; |Scalo 1990 



Rosen & Bregman 1995), implying dumpiness is nec- 



essary to include in any radiative transfer model in- 
volving dust in galaxies. 

In discussing the effects of dust in galaxies, it is 
important to distinguish between a galaxy's attenu- 
ation curve and the dust's extinction curve in that 
galaxy. The dust's extinction curve is directly related 
to the physical properties of the dust grains. An ex- 
tinction curve can be observationally determined by 
using the standard pair method, observing both a red- 
dened and unreddened star of the same spectral type 



and comparing their fluxes (e.g. Massa, Savage, & 



Fitzpatrick 1983; Witt, Bohlin, fc Stecher 1984; Fitz 



patrick 1990). The attenuation curve of a galaxy is 



defined as the change in the galaxy's flux due to the 
presence of dust. This means the attenuation curve 
is dependent on both the physical properties of the 
dust grains and the spatial distribution of the stars 
and dust. Thus, the attenuation curve of a galaxy 
cannot be used directly to probe the composition of 
the dust grains. 

Due to the complexity in studying the dust in the 
integrated SED of a galaxy, multiwavelength obser- 
vations covering a long baseline are needed. For the 
average Milky Way type dust (Ry = 3.1), the dust 
extinction at 1500 A is 2.5 times the extinction at 
5500 A (V) which i s 10 times the extinction a 2.2 jim 
(K) ( |Whittct 1992] ). Therefore, multiwavelength ob- 
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servations probe a large dynamical range of the dust 
attenuation and allow one to better constrain the dust 
spatial distribution and physical properties. 

Multiwavelength observations of starbursts, from 
the UV to the K-band, have recently been used to de- 
rive the general behavior of dust attenuation in star- 
burst galaxies (Calzctti, Kinney, & Storchi-Bergmann 
1994, 1996; |Calzetti 19971 ). Tne average attenuation 
curve for starbursts was determined using a method 
similar to the standard pair method for determining 
extinction curves. The starburst SEDs were collected 
into bins according to their ionized gas E(B — V)i 
values (Storchi-Bergmann, Calzetti, & Kinney 1994; 



Calzctti ct al. 1994) and attenuation curves were de- 
termined using the smallest E(B — V)i bin as the 
unattenuated stellar population. When scaled by 
their respective E(B — V)i values, the attenuation 
curves were roughly similar in shape and strength, 
and this is the behavior of dust in a screen-like geome- 
try. When normalized to an E(B-V) = 1, the resulting 
average starburst attenuation curve has a slope sim- 
ilar to the MW extinction curve, but lacks a 2175 A 
bump, like the SMC extinction curve. Generally, the 
dust distribution can be described by a distribution of 
foreground clumps. The strength of the attenuation 
affecting the stars was found to be different than that 
affecting the gas, implying that the stars have a dust 
covering factor smaller (by ~40%) than the gas. This 
implies the stars and gas have different spatial distri- 
butions. While the flatter attenuation curve can be 
explained as an effect of the radiative transfer in the 
dust, we will argue in this paper the lack of a 2175 A 
bump in the average attenuation curve can only be 
explained by a dust extinction curve with a weak or 
missing bump. 

Without prior knowledge of the dust's type and 
distribution or stellar evolutionary history in star- 
bursts, the only way to investigate the SED of a star- 
burst is to combine a stellar evolutionary synthesis 
model with a model for the radiative transfer in dust 
and simultaneously solve for the characteristics of the 
stars and dust. The aim of this work is to investi- 
gate a large sample of starbursts and derive the ge- 
ometry and type of dust present in starbursts. In 
future work, this will allow us to investigate with 
greater confidence the evolutionary history of indi- 
vidual starbursts. In §0, we discuss our model of a 
starburst galaxy which includes stellar evolutionary 
synthesis and the radiative transfer through the star- 
burst's dust. The observational sample (consisting 



of 30 starburst galaxy SEDs) used in this paper is 
described in §|[ Section ||| presents color-color plots 
useful for detecting dust in starbursts as well as con- 
straining the type and spatial distribution of the dust. 
As an example, we apply our starburst model to the 
galaxies NGC 4385 & NGC 7714 in §§. We discuss 
our results and their implications in §pf Finally, our 
conclusions are presented in §[7| 

2. Starburst Model 

In order to model the SED of starburst galaxies, 
both a stellar evolutionary synthesis model and a dust 
radiative transfer model are needed. We combine 
the stellar evolutionary synth esis mo dels of Leitherer 
& Heckman (1995; hereafter LH95 ) and Bruzual & 
Chariot (1993, 1997; hereafter]BC97]) to obtain the 
intrinsic SEDs of both continuous and instantaneous 
bursts of star formation with ages between 1 x 10 6 
and 15 x 10 9 years. We use a radiative transfer code 
based on Monte Carlo techniques to model the effects 
of dust on the intrinsic SEDs for both MW and SMC 
type dust. 

2.1. Radiative Transfer in Dust 

In order to accurately model the radiative trans- 
fer of the stellar photons through the dust, we rely 



on Monte Carlo techniques (Witt 1977). With these 
techniques, photons are followed through a dust dis- 
tribution and their interaction with the dust is pa- 
rameterized by the dust optical depth (r), albedo, 
and scattering phase function asymmetry (g). The 
optical depth determines where the photon interacts, 
the albedo gives the probability the photon is scat- 
tered from a dust grain, and the scattering phase 
function gives the angle at which the photon scat- 
ters. As the physical properties of dust in different 



environments vary greatly (Witt, Bohlin, & Stecher 


1984 


; |Cardelli, Clayton, & Mathis 1989; 


Fitzpatrick 


1989 


; Mathis & Cardelli 1992 


), we have chosen to 



model two extremely different kinds of dust: MW and 
SMC. For both kinds of dust, we have compiled the 
wavelength dependence (1200-23,000 A) of the opti- 
cal depth, albedo, and g in Table [|. For MW type 
dust, the t/t\/ values arc from Whittct (1992) and 
the albedo and g values were determined from smooth 
curves drawn through many empirical determinations 
taken from the literature (e.g. Gordon et al. 1994 [and 
references therein]; |Calzctti ct al. 1995 ; Lchtincn & 
Mattila 199"^ ). For SMC type dust, the t/t v values 
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are from Prevot et al. (1984) and the albedo and g 
values were kindly computed for us by Sang-Hee Kim 
(Kim, Martin, & Hendry 1994; Kim 1996) assuming 



that the SMC extinction curve is produced by a mix- 
ture of silicate and graphite grains. The characteris- 
tics of the magnitude bands (UV1 to K) are given in 
Table § 

The spatial distribution of the stars and dust was 
investigated by assuming either dusty or shell geome- 
tries. We expect these geometries to bracket the large 
range of observational star/dust distributions. Fol- 
lowing WTC, the dusty geometry has stars and dust 
uniformly mixed through the entire sphere and the 



shell geometry (nuclei geometry in WTC) has a dust- 
free sphere of stars (r < 0.3r max ) surrounded by a 
star-fr ee shell of dust (0.3r max < r < r max ). Unlike 
WTC , we also allow the dust to be either homoge- 
nco us or clumpy (two-phases) on a local scale (but 
see Witt, Gordon, fc Madsen (1997) for an extension 
of WTC] to include clumpy dust). The two-phase 
clumpy dust is characterized by the ratio of inter- 
clump to clump medium density (kz/ki), the filling 
factor (ff) of the high density clumps, and the size 



of the clumps compared to the system size (Witt & 
GoJ-don 199(H). 



The radiative transfer model has the following in- 
puts: the geometry (dusty or shell), homogeneous or 
clumpy dust distribution (fc 2 /fci, ff, & clump size), 
radial r in the V band (equivalent to dust mass), and 
either MW or SMC type dust grain properties (r, 
albedo, and g [see Table [Tj]). The different model runs 
are listed in Table p|. All model run s were computed 
for ff = 0.15 ( |Witt fc Gordon 1996| ), a clump size of 
10% of the system size (radius), and 0.05 < ry < 40. 

2.2. Stellar Evolutionary Synthesis 

The description of the stellar evolution of the star- 



bursts is derived from the work of LH95 and BC97 



LH95 concentrated on the stellar evolutionary syn- 
thesis of massive stars, including nebular emission as- 
sociated with massive stars. They produce SEDs for 
the two different extreme cases of star formation, in- 
stantaneous and constant. For both cases of star for- 
mation, we use an initial mass function (IMF) with 
a Salpeter (1955) slope (a — 2.35), a mass range of 
1-100 Mq, and solar metallicity. We have adjusted 
the absolute level of the SEDs to re flect a mass range 
of 0.1-100 Mq using equation 3 in LH95 . Since their 
model lacks information about the later stages of stel- 



lar evolution, we limit our use of their SEDs to those 
with ages less than 2 x 10 7 years. 

BC97 present a stellar evolutionary synthesis code 
which produces SEDs for both instantaneous and con- 
stant star formation. The SEDs are formed using an 
IMF with a Salpeter (1955) slope (a — 2.35), a mass 
range of 0.1-125 Mq, and solar metallicity. As this 
code lacks information on the nebular continuum, we 
limit our use of these SEDs to the ages between 5 x 10 7 
and 1.5 x 10 10 years. In the continuous star forma- 
tion case, where the nebular contribution from the 
most recent generation of stars is present, we add the 
nebular continuum calculated from a LH95| model of 
constant star formation for an age of 2 x 10 7 years. 



While both LH95 and BC97 provide SEDs com- 



puted for metallicities other than solar, we have cho- 
sen to use only the solar metallicity SEDs for three 
reasons. First, the average metallicity of the ob- 
served starbursts is relatively high, around 0.5 solar 
(see §||). Second, theoretical models of low metallicity 
stellar populations still have difficulties in reproduc- 
ing all the observational constraints (e.g., the ratio 
of blue/red supergiants, see Maeder & Conti [1994]). 
Third, using model SEDs with a higher metallicity 
than actually observed will result in a lower limit 
on the effects of dust. This is due to the fact that 



the lower the metallicity the bluer the SED (Worthey 



1994) and since effects of dust are to redden the SED, 



using solar metallicity SEDs will give a lower limit on 
the dust's effects. In future work, we will use SEDs 
of the observed metallicities when investigating the 
stellar evolutionary history of individual starbursts. 

The upper mass limit on the IMF between the two 
models is different; 100 M for |LH95| and 125 M 



for BC97 







This is not expected to be important as 
the IMF is steep, the amount of mass between 100 
and 125 M is less than 1% of the total mass, and 
its relative contribution to the ionization and heating 
of the interstellar medium is small. 

Using the above stellar evolutionary models, SEDs 
from 1000-23,000 A were generated with ages between 
1 x 10 6 and 15 x 10 9 years. These SEDs were con- 
volved with photometric bandpasses at UV, optical, 
and IR wavelengths to generate the intrinsic colors 
of the starburst model galaxies. In the UV, 7 square 
bandpasses (UV1-UV7) w ere defined with stan dard 
UV zero magnitude fluxes jWesselius et al. 1982| ). For 
the optical and IR, photo metric bandpasses for UB- 
VRIJHK were taken fromlBessell & Brett (1988)1 and 



Bessell (1990). The zero magnitude fluxes for the UB- 
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Table 1 
Dust Physical Properties 



band 


t/tv 


MW 
albedo 


! 


y 


t/tv 


SMC 
albedo 


g 


UV1 


3.11 


0.60 





.75 


5.00 


0.40 


0.53 


UV2 


2.63 


0.67 





.75 


4.36 


0.40 


0.53 


UV3 


2.50 


0.65 





.73 


3.51 


0.58 


0.54 


UV4 


2.78 


0.55 





.72 


3.20 


0.58 


0.51 


UV5 


3.12 


0.46 





.71 


2.90 


0.55 


0.46 


UV6 


2.35 


0.56 





.70 


2.40 


0.56 


0.37 


UV7 


2.00 


0.61 





.69 


2.13 


0.53 


0.35 


u 


1.52 


0.63 





.65 


1.58 


0.46 


0.34 


B 


1.32 


0.61 





.63 


1.35 


0.43 


0.32 


V 


1.00 


0.59 





.61 


1.00 


0.43 


0.29 


R 


0.76 


0.57 





.57 


0.74 


0.41 


0.26 


I 


0.48 


0.55 





.53 


0.52 


0.38 


0.23 


J 


0.28 


0.53 





.47 


0.28 


0.33 


0.21 


H 


0.167 


0.51 





.45 


0.17 


0.30 


0.23 


K 


0.095 


0.50 





.43 


0.11 


0.29 


0.22 



Table 2 
Model Run Inputs 



run 


geometry 


k 2 /ki 


dust type 


1 


shell 


1.0 


MW 


2 


shell 


0.1 


MW 


3 


shell 


0.01 


MW 


4 


shell 


0.001 


MW 


5 


shell 


1.0 


SMC 


6 


shell 


0.1 


SMC 


7 


shell 


0.01 


SMC 


8 


shell 


0.001 


SMC 


9 


dusty 


1.0 


MW 


10 


dusty 


0.1 


MW 


11 


dusty 


0.01 


MW 


12 


dusty 


0.001 


MW 


13 


dusty 


1.0 


SMC 


14 


dusty 


0.1 


SMC 


15 


dusty 


0.01 


SMC 


16 


dusty 


0.001 


SMC 
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VRI bands were calculated by convolving the appro- 
priate bandpass with the observed spectrum of Vega 
([Tug, White, fc Lockwood 1977|). The calibrated 



spectrum of Vega was multiplied by 1.028 before use 
to account for the fact that Vega's V magnitude is 
0.03 ([Hofflcit k Warren 1991J). The zero magnitude 



fluxes for the JHK bands were taken from (Bcssell 
19£(o|). The A oq , A A, and zero magnitude fluxes for 
the adopted bandpasses are presented in Table 



3. Observations 

For this paper, we use a sample of 30 starburst 
galaxies which have SEDs available in the literature. 
For consistency, we provide a brief summary of the 
characteristics of galaxies and the properties of the 
SEDs. We refer the reader to the original papers 



for the observational details ( 


Kinney et al. 1992 ; 


Mc- 


luade, Calzetti, & Kinney 1995; 


Storchi-Bergmann. 


Kinney, & Challis 1995; palzetti 1997 


)■ 



The galaxies are from the IUE Atlas of Kinney et 
al. (1993) and are characterized by the presence of ac- 
tive star formation in their central regions. The star 
formation activity manifests itself as strong UV SEDs 
and intense line emission at optical and IR wave- 
lengths. The galaxies are nearby (median distance 
60 Mpc for H = 50 km s _1 Mpc -1 ) irregulars or spi- 
rals with disturbed morphologies. The metallicity of 
the starbursts' ionized gas is in the range 0.1-2 solar, 
with a median around 0.5 solar. 



Aperture-matched UV and optical SEDs (Kin- 
ic> et al. 1993 ;J\fccjuade I _Ca^^ 1995| 



StoreriT^c^mlmii^Kmnxy are aval 

able for all 30 galaxies and aperture-matched IR pho- 



The availability of spectral information allowed us to 
assess the importance of the emission line contribu- 
tion to the broad band photometry. The IR images 
(Calzetti 1997) were used to obtain the J, H, and K 
photometry with in an area of 10" x 20" in the E-W 
direction. In the IR, the contribution of the emission 
lines to the broad band photometry is negligible, ex- 
cept for the J-band photometry of NGC 4861, where 
the hydrogen recombination line Pa/3 (12818 A) con- 
tributes 9% of the band flux ( palzetti 1997] ). 

As with the SEDs produced with the stellar evo- 
lutionary models, the SEDs of the 30 starbursts were 
convolved with photometric bandpasses at UV and 
optical wavelengths to determine their UV2-I pho- 
tometry. Prior to this, the SEDs were corrected for 
foreground Galactic extinction, shifted to zero red- 
shift, and cleaned of emission lines. The dereddening 
was done using the average Galactic extinction curve 
( |Whittet 19"92| ) and the values of E{B — V)g listed 
in Table 1 of Calzetti, Kinney, & Storchi-Bergmann 
(1994). The SEDs were shifted to zero redshift using 
the values of z listed in the same table. The emission 
lines were removed as the stellar evolutionary syn- 
thesis model SEDs did not include emission lines and 
we wish to compare the observations to these models. 
The IR photometry was transformed from the CIT 
system ( Elias et al. 1982) ) to the homogeneous system 
defined by the equations in Bessel & Brett (1988). 
The resulting magnitudes from the UV, optical, and 
IR are tabulated in Table |. The uncertainties in the 
tabulated magnitudes are 0.10 for UV2-UV7 (Bohlin 



ton|ictry at J, H, and K (Calzetti 1997) is available for 



et al. 1990| ; [Bohlm 1997| ), 0.10-0.20 for U-I (0.10 for 
those galaxies with both UV7 and U data), and 0.09- 
0.14 for J-K (|Kinncy et al. 1993|; |Mcquade, Calzetti 



19 51 the galaxies. The IUE aquired UV portion of 
the SED was corrected to reflect the newest IUE ab- 



& Kinney 1995; Storchi-Bergmann, Kinney, & Chal- 
lis 1995|j" Calzetti 1997). Due to mismatches between 



solute calibration ( Bohlin 1996 ) . The match in obser- 
vational apertures ensured that the same region was 
sampled within each galaxy. Since the UV SEDs were 
acquired with IUE, apertures of 10" x 20" were used 
for the optical and IR observations. Differences in ori- 
entation between the IUE and optical apertures pro- 
duced an estimated 10-20% flux mismatch between 
the two wavelength regions. This size aperture sub- 
tends large regions within each galaxy, a diameter of 
^4.5 kpc at the median distance of 60 Mpc. 

In most cases, the UV and optical SEDs cover the 
wavelength range 1250-7500 A continuously with a 
resolution 6 A in the UV and ~10 A in the optical. 



the apertures used in obtaining the optical and IR 
photometry, the relative uncertainties between these 
two wavelength regions are 0.10-0.20 magnitudes. In 
all the plots in this paper, we have assumed a un- 
certainty of 0.10 for UV2-UV7 data, 0.10 for U-I for 
galaxies with UV7 and U data, 0.15 for U-I for galax- 
ies without both UV7 and U, and 0.125 for J-K. For 
colors involving only J-K, we have used a uncertainty 
of 0.07 for (J-H) and 0.05 for (H-K) ( palzetti 1997| ). 



4. Color-Color Plots 

The main problem with accurately modeling the 
ages of the stellar populations and effects of dust in a 
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Table 3 
Magnitude Band Definitions 



band 


Acq 

[A] 


AA 

[A] 


F A a 

[ergs cm~ 2 s -1 A -1 ] 


UV1 


1250 


250 


3.63xKT 9 


UV2 


1515 


280 


3.63xl0~ 9 


UV3 


1775 


260 


3.63xl0~ 9 


UV4 


1995 


200 


3.63xl0~ 9 


UV5 


2215 


260 


3.63xl0~ 9 


UV6 


2480 


290 


3.63xl0~ 9 


UV7 


2895 


560 


3.63xl0~ 9 


U 


3605 


640 


3.41xl0~ 9 


B 


4413 


959 


6.60X10" 9 


V 


5512 


893 


3.70xl0~ 9 


R 


6594 


1591 


2.35xl0~ 9 


I 


8059 


1495 


1.15X10" 9 


J 


12369 


2034 


3.12xl0- 10 


H 


16464 


2862 


1.14xl0- 10 


K 


21578 


2673 


3.94xl0" n 



& F\ is the flux corresponding to a magnitude 
of zero. See text for details. 
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Table 4 

Starburst Galaxy Photometry 



name 


UV2 


UV3 


UV4 


UV5 


UV6 


UV7 


U 


B 


V 


R 


I 


J 


H 


K 


IC 1586 


14.27 


14.34 


















14 


77 


15 


32 


14 


92 








13 


58 


12 


83 


12.52 


HARO 15 


12.99 


13.14 


















14 


37 


15 


02 


14 


78 


14 


53 


14.05 


13 


54 


12 


92 


12.58 


MRK 357 


13.29 


13.55 


















14 


90 


15 


87 


15 


79 








14 


55 


13 


86 


13.39 


IC 214 


14.10 


14.14 


















15 


59 


16 


03 


15 


48 








12 


66 


11 


85 


11.37 


NGC 1140 


11.96 


12.17 


















13 


11 


13 


75 


13 


43 


13 


18 


12.74 


12 


54 


11 


86 


11.62 


NGC 1510 


12.91 


13.14 


13 


42 


13 


56 


13 


87 


14 


02 


14 


33 


14 


87 


14 


55 


14 


25 














NGC 1569 


9.68 


9.72 


9. 


95 


10 


18 


10 


19 


10 


46 


10 


71 


11 


66 


11 


47 


11 


36 




10 


93 


10 


29 


10.12 


NGC 1614 


14.16 


14.16 


















14 


18 


14 


57 


13 


91 


13 


36 


12.57 


11 


46 


10 


60 


10.07 


NGC 1705 


10.55 


10.94 


11 


23 


11 


40 


11 


82 


12 


27 


12 


70 


13 


61 


13 


43 


13 


22 


12.79 












NGC 3049 


13.54 


13.58 


13 


80 


13 


92 


14 


08 


14 
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Fig. 1.— The color-color plane (U-B) vs. (UV2-V) is 
plotted for both the starburst observations and mod- 
els. The dust reddening trajectories can start from 
any stellar evolutionary synthesis model, but to plot 
them, they were attached to an arbitrary stellar evolu- 
tionary synthesis model. The clumpy geometry refers 
to a kijk\ = 0.01. The words burst and const in the 
upper legend refer to the kind of star formation (in- 
stantaneous or constant, respectively) assumed in the 
stellar evolutionary synthesis models. The age vector 
points in the direction of increasing age of the stellar 
evolutionary synthesis SEDs. 



galaxy is that aging the stellar population and adding 
dust both redden the galaxy's SED. A non-exhaustive 
search of the many color-color plots turned up a hand- 
ful of plots where the effects of stellar population 
changes were in a different direction in the color- 
color plane than those due to the dust. By exam- 
ining where the sample of observed starburst galax- 
ies fall on these color-color plots, information about 
the dust (type/geometry) and the stellar population 
(burst/constant and age) can be derived. 

The color-color plot (U-B) vs. (UV2-V) is plotted 
in Figure |l| and plainly shows the presence of dust 
in starburst galaxies. All of the stellar evolutionary 
synthesis models fall in a well defined line, with age in- 
creasing downward, while the actual observations all 
lie to the right of this line. Combining different stel- 
lar evolutionary models to approximate a starburst 
with a complex stellar evolutionary history will not 
move the starburst off of the line defined by the sim- 
ple stellar evolutionary models. All of the reddening 
trajectories point to the right but some are not long 



enough to explain the spread in the observed star- 
bursts. Those which cannot explain the spread are 
the dusty geometry models which assume the same 
global distributions of stars and dust, be they homo- 
geneous or clumpy. The shell geometry can explain 
the observed spread and we adopt this as the geome- 
try for starburst galaxies. This geometry results in a 
screen-like behavior of the dust which allowed Calzetti 
et al. (1994) to derive an average attenuation curve for 
starbursts in a manner analogous to the standard pair 
method. It is important to remember that an attenu- 
ation curve is an extinction curve convolved with the 
dust geometry and thus an attenuation curve cannot 
be used to probe the type of dust grains directly. 

From Figure ^, the stellar population can be con- 
strained by moving along the reddening trajectories 
back to the line defined by the stellar evolutionary 
models. The stellar populations seem to arise from 
burst star formation models with ages between 8 x 10 6 
and 2 x 10 8 years or constant star formation models 
between 1 x 10 8 and 15 x 10 9 years. Actually de- 
termining the star formation history of an individual 
galaxy is quite complex and requires fitting the SED 
from the UV to the IR. See §|| for an example of how 
this can be done. The absence of starbursts (from our 
UV selected sample) with ages < 8 x 10 6 years sug- 
gests the initial phase of starbursts is hidden within 
dust clouds and is difficult to observe in the UV. In 
addition, the age we derive is the average age of the 
starburst regions falling in the large aperture used to 
obtain the SEDs. Thus, starburst regions with ages 
< 8 x 10 6 years could be present, but washed out by 
the presence of older starburst regions. 

Figures show examples of useful color-color 
plots (one each for the UV, optical, and IR) which 
show differences between stellar population changes 
and dust reddening. The (UV2-UV3) vs. (UV2-UV7) 
color-color plane (Figure ^) is a good plane in which 
to investigate whether SMC or MW type dust best 
matches the type of dust in starbursts. The observa- 
tions fall between the SMC and MW reddening tra- 
jectories, but closest to the SMC reddening trajecto- 
ries. This implies the far-UV slope of the starbursts' 
extinction curve has a value bet ween the MW and 
SMC, but closest to the SMC. In §0, we investigate 



the starbursts' UV extinction curve in more detail, 
focusing on the 2175 A bump. The optical (U-B) 
vs. (B-V) color-color plane (Figure ||) does show that 
dust is present, but the effects are not as definitive as 
those shown in the previous two figures. This points 
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Fig. 2.— The color-color plane (UV2-UV3) vs. (UV2- 
UV7) is plotted. The number in parenthesis after the 
word clumpy in the lower legend refers to the value 
of k2/k%. Otherwise, this figure was constructed in a 
manner similar to Figure 1. 



to the need, in investigating the dust in starbursts, for 
SEDs covering more than just the optical. The (J-H) 
vs. (H-K) color-color plane (Figure ^) also shows the 
effects of dust, but some degeneracy exists between 
changes in age and the effects of dust. In addition, 
not all of the starbursts can be traced back to a stellar 
population model through the reddening trajectories. 
This is probably due to problems with stellar evolu- 
tionary models in the IR region (see LH95 and Lancon 
& Rocca-Volmerange [1996] for details). While all of 
the above color-color plots (Figures |l]-[|) clearly show 
the presence of dust in starbursts, the diagnostic util- 
ity of color-color plots is greatest in the UV. 



4.1. Absence of the 2175 A bump 

Due to its prominence, the 2175 A extinction bump 
is an excellent probe of the type of dust in a galaxy. 
This is illustrated by the average extinction curves of 
the MW, LMC, SMC, and M31, all of wh i ch show 
different 2175 A strengths ( Fitzpatrick 1989] ; Bianchi 
t al. 1996| ). The lack of a depression at 2175 A in 



SEDs of starbursts has raised the question: does the 
dust in starbursts lack a 2175 A bump or are there 
radiative transfer effects associated with the spatial 
distribution of the dust which weaken the bump be- 
low observational detection? To investigate this ques- 
tion, two color-color plots were constructed. Figures ^ 
& | show the (UV4-UV5) vs. (UV4-UV6) and (U-B) 




Fig. 3. — The color-color plane (U-B) vs. (B-V) is 
plotted. The number in parenthesis after the word 
clumpy in the right legend refers to the value of kifk\. 
Otherwise, this figure was constructed in a manner 
similar to Figure 1. 
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Fig. 4. — The color-color plane (J-H) vs. (H-K) is 
plotted. The number in parenthesis after the word 
clumpy in the lower legend refers to the value of /k\ . 
Otherwise, this figure was constructed in a manner 
similar to Figure 1. 
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Fig. 5.— The color-color plane (UV4-UV5) vs. (UV4- 
UV6) is plotted. The number in parenthesis after the 
word clumpy in the lower legend refers to the value 
of kz/kx. Otherwise, this figure was constructed in a 
manner similar to Figure 1. 




-0.6 -0.4 -0.2 0.0 0.2 

UV4 - UV5 



Fig. 6.— The color-color plane (U-B) vs. (UV4-UV5) 
is plotted. The number in parenthesis after the word 
clumpy in the lower legend refers to the value of kifkx- 
Otherwise, this figure was constructed in a manner 
similar to Figure 1. 



vs. (UV4-UV5) color planes, respectively. The place- 
ment of the three bandpasses of UV4, UV5, & UV6 
(A cq = 1995 A, 2215 A, & 2480 A) make them ideal 
for studying the 2175 A bump. The color (UV4-UV6) 
measures the slope underneath the 2175 A bump and 
the color (UV4-UV5) is sensitive to the short wave- 
length slope of the 2175 A bump. 

Both figures clearly show the SMC reddening tra- 
jectories generally point toward the starburst obser- 
vations, while the MW reddening trajectories com- 



dust in the SMC than in the MW. This agrees with 
the notion that the dust in starbursts would be simi- 
lar to the weak 2175 A bu mp 30 Dor regi on dust, the 



mini-starburst next door (Walborn 1991) 



In Figure M, there are a number of starbursts which 
cannot be traced back to the stellar evolutionary mod- 
els through a reddening trajectory. This is not sur- 
prising as we have only investigated two extinction 
curves and the extinction curves of dust in the MW, 
LMC, and SMC are known to possess large variations 



pletely miss the observations. In Figure [sj the MW ( 
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from the observed starbursts. One of the MW red- 
dening trajectories shows how radiative transfer ef- 
fects can reduce the strength of the 2175 A bump. 
In Figure ||, the reddening trajectory for model run 
4 (shell, MW clumpy with k 2 /kx = 0.001) loops back 
toward the stellar evolutionary models. This is caused 
by the output spectrum being dominated by unatten- 
uated stars which suppress both the 2175 A bump 
as well as the slope of the attenuation curve in the 
UV. While the starbursts do not show the presence 
of the 2175 A bump, they do show the signature of 
a substantial UV slope in their attenuation curves. 
Therefore, radiative transfer effects alone cannot ex- 
plain the lack of a 2175 A bump in starbursts and, 
so, the dust in the average starburst is more like the 



bursts, we would need an extinction curve in which 
the color excess in (UV4-UV5) is greater for a given 
color excess in (UV4-UV6) than that in the SMC 
extinction curve. This implies an extinction curve 
with a curvature between UV4 and UV6 greater than 
curvature in the same region in the SMC extinction 
curve. Confirmation of this speculation awaits future 
work where we plan to individually fit each of the 
starbursts and derive the average starburst extinction 
curve. 

5. NGC 4385 & NGC 7714 

Fitting an individual galaxy's SED is difficult as 
there are usually many starburst models (combina- 
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tions of stellar evolutionary synthesis and radiative 
transfer in dust) which can produce the same SED. As 
an example of this degeneracy, we fit the two galaxies 
NGC 4385 and NGC 7714 for which we have com- 
plete SEDs from UV2 to K. To test the goodness of 
the fit for a specific starburst model, we used the x 2 
(reduced x 2 ) criterion after first normalizing the ob- 
served and model SEDs at the K band. The x 2 was 
computed using 



x 2 



1 

i=l 



F {Xj) - F m (Xj) 



(1) 



where n is the number of points in the observed SED, 
d is the number of degrees of freedom, F (Xi) is the 
observed flux at the ith point in the SED, F m (Xi) 
is the same for the model flux, and_o(Aj]_Js the ith 
uncertainty in the observed flux ( Taylor 1982 ). The 
number of degrees of freedom d = (n — 1) as both 
the observed and model SEDs were normalized at 
the K band. Starburst models with x 2 > 1.8 were 
rejected as the probability of such models agreeing 
with the observations was small. With this limiting 
value of x 2 1 we wm collect 96.3% of the starburst 
models which fit the observed SED within the uncer- 



tainties (Taylor 1982). Such fits were obtained for 17 
and 46 different starburst models for NGC 4385 and 
NGC 7714, respectively. 

The parameters of the starburst models which fit 
are interesting. Only models containing dust with a 
SMC extinction curve fit the data, consistent with 
Using our results from §[| and 



our findings in §4.1 



restricting the fits to models with a shell geometry, 
the number of starburst models which fit are 8 for 
NGC 4385 and 18 for NGC 7714. In order to choose 
among the starburst models which fit the observed 
SEDs, more information about each starburst (e.g. 
the number of ionizing photons) must be used and 
we plan to do this in a future paper. The best fit 
starburst model for NG C 4385 (x 2 = 1.11) was for a 
constant star formation BC97 model with an age of 
7 x 10 9 years, a shell geometry with ry = 0.75, and 
clumpy (fc 2 /fci = 0.001) SMC type dust. This best 
fit starburst model, along with the observed SED of 
NGC 4385 is displayed in Figure @a. The best fit 
starburst model for NG C 7714 (x 2 = 0.49) was for a 
constant star formation BC97 model with an age of 
1 x 10 9 years, a shell geometry with Ty = 1.00, and 
clumpy (k 2 /ki = 0.001) SMC type dust. This best 
fit starburst model, along with the observed SED of 
NGC 7714 is displayed in Figure 0b. 
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Fig. 7.— The best fit starburst model for NGC 4385 
(a, left) and NGC 7714 (b, right) are plotted along 
with their observed SEDs. 
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To illustrate the difference between extinction and 
attenuation curves, Figure || displays the SMC ex- 
tinction curve and the average starburst attenua- 
tion curve (Calzetti et al. 1994) with the attenua- 
tion curves from the best fit starburst models for 
NGC 4385 and NGC 7714. For both galaxies, the 
attenuation curves are flatter than the dust extinc- 
tion curve and this is consistent with the Calzetti et 
al. (1994) attenuation curve. Also, the shape of the 
attenuation curve is different between NGC 4385 & 
NGC 7714 even though both curves are for a shell ge- 
ometry with a clumpy {k^l^-i = 0.001) distribution of 
SMC dust. This neatly shows the dependence of the 
shape of an attenuation curve on both the geometry 
and amount of dust. The systematic difference be- 
tween the shape of both the NGC 4385 and NGC 7714 
attenuation curves and the Calzetti et al. (1994) at- 
tenuation curve is likely due to differences between 
the actual dust albedo and the adopted albedo (Ta- 
ble@). 

Calzetti et al. (1994) determined the effective B 
band optical depth (r B , c ff) to be 0.62 for NGC 4385 
and 0.42 for NGC 7714. These values of t b . cS were 
determined from observations of Ha and H/3; there- 
fore they measure the TB iG ff of the starburst's gas. 
The starburst models in this paper measure the star- 
burst's stars and from the best fit models the TB, e ff 
was 0.26 for NGC 4385 and 0.33 for NGC 7714. This 
result agrees with the finding by Calzetti (1997) that 
the gas in starbursts shows a larger dust attenuation 
than the stars. The gas is likely associated with O 
stars (H II regions), while the stellar continuum has 
contributions from a wide range of stars. As O stars 
are short lived, on average they will be more deeply 
embedded in their birth clouds and thus the gas will 
show a larger dust attenuation than the stellar con- 
tinuum. 

6. Discussion 

The discovery that the dust in starbursts is closer 
to SMC type dust than MW type dust is not totally 
unexpected ( palzetti et al. 1994 ; Mas- Hesse & Kunth 
1996). It has been known that the attenuation curve 
of starbursts is flat and does not show the presence 
of the 2175 A bump, but as a result of this work it 
is clear this is due to SMC-like dust grains and not 
radiative transfer effects. The 30 Dor region of the 
LMC forshadows this conclusion as the dust associ- 
ated with the massive star formation in the 30 Dor 
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Fig. 8. — This plot shows a comparison of the SMC 
extinction curve and the Calzetti et al. (1994) aver- 
age starburst attenuation curve with the attenuations 
curves from the best fit starburst models of NGC 4385 
and NGC 7714. The average starburst attenuation 
curve was created by normalizing the Calzetti et al. 
(1994) starburst attenuation curve to the stellar con- 
tinuum (E(B-V) = 1) and adding R v = 3.1. The 
curves are normalized to the V band. The effective 
Ty for NGC 4385 was 0.19 and the input Ty was 0.75. 
The effective t v for NGC 7714 was 0.25 and the in- 
put Ty was 1.00. The effective Ty is defined as the Ty 
an overlaying screen would need to produce the same 
reduction in flux. 
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region is very different than the dust associated with 
the rest of the LMC. The 30 Dor region dust has a 
weaker 2175 A bump and stronger far-UV rise than 
the dust in the rest of the LMC. Like the 30 Dor re- 
gion, it is unlikely the dust in starbursts formed with 
a weak or negligible 2175 A bump and strong far-UV 
rise. It is more likely that the dust in starbursts has 
been modified by the starburst activity to more SMC- 
like in nature. Again, the 30 Dor region has shown 
i ts identity as t he Rosctta Stone for starburst galaxies 
dWalborn J99j| ). 

With this conclusion in hand the question arises: 
How does starburst activity modify dust such that 
it loses its 2175 A bump and strengthens its far-UV 
rise? The answer lies in the different grain popula- 
tions which taken together explain dust extinction 
curves. Dust grains are thought to be composed of 
three different populations: classically sized grains 
which produce the optical- IR extinction, small grains 
which produce the 2175 A bump, and small grains 



which produce the far-UV extinction rise ( Mathis 
199 5 |Zubko, Krclowski, We gner 199q; |Li fe Green 



ber i 1997). The SMC-likc extinction curve in star- 



bursts implies the 2175 A bump grains are destroyed 
and the number of far-UV rise grains have been in- 
creased. The 30 Dor region has less intense star for- 
mation leading to a similar but a notably different 
extinction curve. This extinction curve implies that 
the 2175 A bump grains are still present, but in re- 
duced numbers, and the far-UV rise grains have in- 
creased their numbers. Finally, the Orion region has 
a low level of star formation (~1% that of 30 Dor 
Hunter et al. 1995]) which has also affected the dust's 
extinction curve. The Orion region extinction curve 
has a reduced 2175 A bump and a nearly fiat far-UV 



rise (Cardelli & Clayton 1988). This implies both the 
2175 A bump and far-UV rise grains have been re- 
duced in number. 

There seems to be a trend in the dust extinction 
curve with star formation activity. This can be ex- 
plained by the relative importance of shocks and stel- 
lar radiation in the modification of dust grains. In 
the Orion region, the stellar radiation pressure and 
evaporation of du st grains hav e selectively removed 
the small grains ( McCall 1981 ; Cardelli fc Clayton 
19^||). In starbursts, the dust modification is most 
likely dominated by supernovae shocks which are ef- 
ficient at shattering large grains, thus increasing the 
small grain population ( Jones, Ticlcns, fc Hollenbach 
199q). This explains the increase in the number of 



far-UV rise grains, but the absence of the 2175 A 
bump grains implies they do not survive the shocks 
and are either vaporized or shattered. The 30 Dor 
region lies between the Orion region and starburst 
galaxies in star formation activity and, thus, has an 
extinction curve intermediate between the two. The 
trend with star formation activity implies the mate- 
rial which makes up the 2175 A bump grains is not 
as robust as the material which makes up the far-UV 
rise grains. The identity of material which gives rise 
to the far-UV rise and, especially, the 2175 A bump 
has remained elusive for many years. For this reason, 
we do not attempt to discuss the material which pro- 
duces the 2175 A bump, but only to note that the 
trend with star formation activity of the far-UV rise 
and 2175 A dust grains is a clue to the identity of the 
material making those dust grains. 

The similarity of the average extinction curves of 
starbursts and the SMC deserves comment. The lack 
of a 2175 A bump and the steep far-UV rise in the 
SMC extinction curve is attributed to the SMC's low 



meta ll icity ( Prevot et al. 1984 ; Clayton fc Martin 



1985 ; Fitzpatrick 1986 ). Since our sample of star- 
bursts possess a large range of metallicities, 0.1-2 so- 
lar (see section ||), this is unlikely to explain the SMC- 
like extinction curve we found in starbursts. We sug- 
gest the SMC-like extinction curve in starbursts is 
a result of the massive star formation. Thus, both 
metallicity and star formation activity may influence 
dust properties significantly. 

Observations which resolve starburst galaxies show 
that the star formation occurs in localized regions 
analogous to the star clusters and associations of 



our galaxy (p'Connell, Gallagher, & Hunter 1994 


O'Connell et al. 1995 




Whitmore & Schweizer 1995 


Vacca 1996 


). Not all starbursts seem to form stars 



a global starburst. These observations allow us to 
interpret the fact that we were able to model the ef- 
fects of dust in starbursts with a radiative transfer 
model with a clumpy shell geometry. The Ty deter- 
mined from our starburst model is interpreted as the 
flux weighted average Ty to many individual star clus- 
ters. From this, it is not surprising we find fairly low 
Ty's as star clusters possessing a large Ty would not 
contribute much to the starburst SED. The shell ge- 
ometry is a realistic approximation of the geometry 
of the individual star clusters. Near the center of the 
star clusters where the radiation field is intense, dust 
grains would be unable to survive and, thus, the inner 
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region would be fairly free of dust. 

We were able to fit the SEDs of NGC 4385 & 
NGC 7714 by using the same radiative transfer model 
and only varying the amount of dust. This particular 
model run (#8, Table |) had SMC dust, a shell geom- 
etry, a clumpy dust distribution with kz/ki = 0.001. 
This gives hope that it will be possible to investi- 
gate the stellar evolutionary history of individual star- 
bursts as well as the type and amount of dust present. 
These fits were done using only the broad band SED 
of the starbursts and a proper fit would use the full 
SED. We plan to fit the full SED and fit both the 
continuum and spectral lines in future work where we 
plan to investigate each starburst individually. 

7. Conclusions 

From the comparison of 30 starburst SEDs with a 
model for starbursts which included both stars and 
effects of dust, we concluded the following: 

1. The dust in starbursts has an extinction curve 
similar to that found in the SMC, i.e. lacking a 
2175 A bump and having a far-UV rise with a 
slope intermediate between the MW and SMC. 
The dust is unlikely to have been formed with 
these properties, but to have acquired them as 
a result of the massive star formation occurring 
in the starburst. 

2. The spatial distribution of stars and dust is 
best described by a shell geometry. This ge- 
ometry has an inner dust-free sphere of stars 
surrounded by a star- free shell of dust. 

3. The dust distribution in starbursts is best de- 
scribed by a locally clumpy distribution. 

4. The age of the stellar populations in this UV 
selected sample of starbursts ranges between 8 x 
10 6 and 2 x 10 8 years if a burst star formation 
is assumed and 1 x 10 8 and 15 x 10 9 years if 
constant star formation is assumed. 

5. The effects of dust on the SED of a starburst are 
characterized by an attenuation curve which is 
SMC-like extinction curve convolved with the 
radiative transfer effects associated with a shell 
star/dust geometry with a clumpy dust distri- 
bution. The wanton use of a SMC extinction 
curve to correct a starburst 's SED is not sup- 
ported by this work. 
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